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ABSTRACT 

o 

' The wavelength dependence of interstellar polarization due to dust in M31 has been 

' observed along four sightlines. Only one sightline had been measured previously. The 

■ globular clusters, S78, S150, S233 and Baade 327 were used as point sources to probe 

the interstellar dust in M31. The Serkowski law produces good fits for all the sightlines 
although the relationship between K and Xmax niay be different from that found in the 
Galaxy. The results of this study imply that the slope K/Xmax may be significantly 
larger in M31. The Serkowski curves are significantly narrower than those of the same 
' Xmax in the Galaxy and may require extreme modifications to the size distributions of 

^ ■ silicate particles. The fits for the four sightlines reveal values of Xuax ranging from 4800 

(N ■ to 5500 A. These are consistent with average values of Xmax measured in the Galaxy 



o 



and the Magellanic Clouds. The range measured for M31 corresponds to Ry values of 



Tj" . 2.7 to 3.1. The range in Ry seen in the Galaxy is 2.5 to 5.5 implying, for this small 

sample, that the average size of interstellar grains in M31 is typically smaller than that 



seen for Galactic grains if the nature of the grains is the same. Also, the polarization 
efficiency for these sightlines is large although some bias is expected since sightlines 



O 

$-H ' known to have significant interstellar polarization were selected for the sample. 

a 

Subject headings: M31, interstellar dust, polarization 

■ 

•i-H - 

5t , 1. Introduction 

There is an average Milky Way extinction relation, depending on only one parameter, the ratio 
of total-to-selective extinction, Ry, which is applicable to a wide range of interstellar dust environ- 
ments (Cardelli, Clayton, & Mathis 1989 (CCM)). Similarly, there is an interstellar polarization 



^ A portion of the results presented here made use of the Multiple Mirror Telescope Observatory, a facility operated 
jointly by the University of Arizona and the Smithsonian Institution. 

^Department of Physics & Astronomy, Louisiana State University Baton Rouge, LA 70803; email: gclay- 
t on@ fenway . phys.lsu.edu 

^Space Science Institute, 1540 30th Street, Suite 23, Boulder, CO 80303-1012; email: wolff@colorado.edu 

Steward Observatory, University of Arizona, Tucson, AZ 85721; E-mail: kgordon, psmith@as.arizona.edu 
^Space Astronomy laboratory, University of Wisconsin, Madison, WI 53706; email: khn, brian@sal.wisc.edu 



- 2 - 



relation, correlated with a single parameter, Xmax, the wavelength at which the interstellar polar- 
ization is a maximum (Serkowski, Mathewson, & Ford 1975; Clayton et al. 1995; Martin, Clayton, 
& Wolff 1999). The parameters, Ry and Xmax, both involve integrals over the size distribution 
of interstellar dust grains and are well correlated with each other (Serkowski et al. 1975; Whittet 
&; van Breda 1978; Clayton & Mathis 1988). The existence of these relations, valid over a large 
wavelength interval, suggests that the environmental processes that modify the grains are efficient 
and affect all grains. However, the CCM relation is not a universal law. It does not fit the observed 
extinction in the sightlines observed in the Magellanic Clouds (e.g., Clayton et al. 1996; Gordon & 
Clayton 1998; Misselt, Clayton, &: Gordon 1999; Gordon et al. 2003). The differences in observed 
extinction properties seen among the Local Group galaxies may be due to various environmen- 
tal factors including metallicity. The metallicity of the LMC and SMC is lower than the Galaxy. 
Therefore, the nature of the dust in M31 is of great interest because its metallicity is similar to 
that in the Milky Way. 

The nature of interstellar dust in M31 is summarized in Lequeux (2000). The UV extinction 
data for M31, which are based on only a few lightly reddened sightlines, show a CCM-Iikc far-UV 
extinction and a weak 2175 A bump (Bianchi et al. 1996). Several studies have suggested that the 
extinction in M31 in the U-band differs from the Galaxy (Searle 1983; lye &; Richter 1985; Massey 
et al. 1995). It has also been suggested that there is a systematic change in the U-band extinction 
with distance from the center of M31 (lye & Richter 1985; Preedman &; Madore 1990). Various 
extimates of Ry ranging from 2.5 to 3.8 have been made using individual stars and clusters in M31 
(Kron & Mayall 1960; Sharov & Lyutyj 1989; Nedialkov & Vcltchcv 1999). Recently, Barmby et 
al. (2000) find that the U- to K-band extinction curve in M31, derived from a large number of 
globular clusters, is consistent with CCM with Rv'=3.1. 

The interstellar polarization in the LMC and SMC shows a similar wavelength dependence 
to that seen in the Milky Way and can be well fit by the Serkowski relation (Clayton, Martin, 
&; Thompson 1983; Clayton et al. 1996; Rodrigues et al. 1997). Some interstellar polarization 
measurements in M31 have been made using globular clusters as point sources (Hiltner 1958; 
Martin & Shawl 1979, 1982). However, the wavelength dependence of interstellar polarization has 
been measured along only one line of sight (Martin & Shawl 1979). This polarization measurement is 
consistent with the Serkowski relationship and a small value of Xmax- We report further observations 
of the wavelength dependence of interstellar polarization in M31. 

2. Observations 

Optical spectropolarimctry was obtained using the 6.5 m telescope of the Multiple Mirror 
Telescope (MMT) Observatory located on Mount Hopkins, Arizona. Observations were made of 
four globular clusters listed in Table 1 on 2003 October 30 UT. All of the observations made use 
of the CCD Spectropolarimeter (Schmidt, Stockman, &; Smith 1992), upgraded with a 1200 x 800 
pixel, thinned, antireflection-coated, UV-sensitized CCD, and an improved camera lens and half- 
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wave plate. All observations utilized a 600 line mm ^ grating for high throughput and wide spectral 
coverage, 4100-8300 A. The FWHM resolution with this grating and a slit width of 1" is 17 A 
(3 pixels). A polarimetric measurement sequence involves four separate exposures that sample 16 
orientations of the semiachromatic half-wave plate and totalled 480-1440 s of integration for each 
of the four clusters through light cloud cover. Polarization position angles were referenced to the 
equatorial system using interstellar polarization standards (Schmidt, Elston &; Lupie 1992). The 
spectral flux distribution is also obtained in the course of the observations and is calibrated relative 
to spectrophotometric standard stars selected from the IRAF database (Massey et al. 1988). The 
data have been binned and are plotted in Figure 1 along with the spectra of the globular clusters. 

In addition, spectropolarimetric observations of S78 were obtained with the HPOL spectropo- 
larimeter on the 3.5 m Wisconsin-Indiana- Yale-NOAO (WIYN) telescope on Kitt Peak. The spec- 
tropolarimeter is a modified Boiler & Chivens spectrograph with a rotating superachromatic wave- 
plate after the slit (Nordsieck & Harris 1996; Wolff, Nordsieck, & Nook 1996). Spectra are obtained 
at eight waveplate angles on a Reticon 400 x 1200 CCD. Two gratings are used, producing spectra 
and 3 pixel resolutions of 3200-6200 A (9 A) and 6000-10500 A (12 A), respectively Two polarized 
star spectra and two polarized sky spectra are obtained simultaneously. The calibration accuracy 
in polarization is ibp/500 and in position angle is ±1°. There is an uncertainty of 0.1% in the cor- 
rection for instrumental polarization. This small uncertainty does not affect the results presented 
here. The WIYN spectropolarimetric data were obtained in 1996 December and 1998 October. The 
four nights of data on S78 have been combined and binned to increase the signal-to-noise. They 
are plotted in Figure 1. 

The polarization measured toward M31 must be corrected for a foreground component. These 
lines of sight pass through the disk and halo of the Galaxy, and then the halo and disk of M31. 
However, since M31 lies at a Galactic latitude of -22°, the foreground extinction due to dust in 
the Galaxy is small. Several estimates have been made of the reddening foreground to M31. They 
all indicate that E(B-V) ~ 0.08-0.11 (Bianchi et al. 1996 and references therein). Using the 
well-known Galactic relationship between reddening and interstellar polarization, the maximum 
foreground polarization is ~ 0.9% (Serkowski et al. 1975). Martin &: Shawl (1979) estimated the 
actual foreground polarization by measuring a sample of Galactic stars lying in the direction of 
M31 as well as the nucleus of M31 itself. The foreground dust is patchy so there is no way to know 
a priori what the foreground component is along any given sightline. So, Martin & Shawl (1982) 
adopted a best-guess value of Pma3;=0.4%, Xmax= 5500 A with a position angle of 90° by averaging 
several foreground sightlines. 

The new observations reported here have been corrected for foreground polarization in the 
same manner as Martin k, Shawl (1982). A Galactic Serkowski curve, as described above, was 
assumed and the foreground polarization at each wavelength was calculated and subtracted. The 
foreground-corrected polarizations are plotted in Figure 1. The observed and foreground corrected 
values are listed in Tables 2-6. The new data are consistent with the Martin &: Shawl (1982) data 
for S78 in the overlapping wavelength range. Similarly, the three wavelength bins for S150 and one 
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bin for S233, agree with our new data within 2a. There is only one filterless observation of Baade 
327 having low S/N. The Martin k. Shawl filter data are plotted for comparison in Figure 1. 

3. Discussion 

The interstellar polarization due to dust in M31 was first measured by Hiltner (1958) using 
globular clusters as point sources. He made unfiltcrcd observations of about 20 clusters finding 
several with significant amounts of polarization. Filter polarimetry of a similar number of clusters 
was obtained by Martin &; Shawl (1979,1982). Only one sightline in their sample, toward S78, 
had enough wavelength coverage for a Serkowski fit to be made. That sightline showed Pmaa; = 
1.98±0.07 %, and Xmax = 3500±300 A using the Wilking et al. (1980) formulation of the Serkowski 
relation. 

Our new data for S78 and three additional clusters cover approximately the same wavelength 
range (4100-8300 A) as Martin k, Shawl but with many more data bins and smaller uncertainties. 
We performed nonlinear least-squares fits to the foreground corrected data for all four clusters 
using the Serkowski law (Serkowski et al. 1975; Wilking et al. 1980; Wilking, Lebofsky, & Rieke 
1982; Whittet et al. 1992). We found that the K= l.^^Xuax relationship, that holds for Galactic 
interstellar polarization, does not produce good fits for the M31 sightlines (Whittet et al. 1992). 
However, the Serkowski law with K as a free parameter does produce good fits for all of the 
sightlines. The three fit parameters are listed in Table 1. The relationship between K and Xmax is 
shown in Figure 2. Although, we do not have enough sightlines to define the ¥^-Xmax relationship 
for M31, Figure 2 shows that the slope may be steeper than in the Galaxy. 

The Y^-Xmax relationship for the M31 globulars is very interesting in that it produces distinctly 
narrower "Serkowski curves" than those of the same Xmax in the Galaxy. As has been shown many 
times in the past (Wolff et al. 1994; Clayton et al. 1995; Martin et al. 1999), grain models have 
generally been able to reproduce observed Galactic interstellar polarization curves by small changes 
to the size distribution of aligned grains, alignment efficiency, upper integration limit, etc. This is 
definitely not the case for the M31 curves in that similar changes to a distribution of astronomical 
silicates failed to produce K values of the observed amplitude. Even more extreme modifications 
to the size distribution of silicate particles (use of "flat" distributions, very narrow size ranges, 
etc.) tended to produce K values that still remained 30-40% too low, at best. As a result, it would 
appear that there is something intrinsically different (and interesting) about the grains producing 
the observed polarization in M31 from those in our Galaxy. Wc intend to investigate this behavior 
in much greater depth after our sample of M31 globulars has been increased. 

Previously, only S78 had enough wavelength coverage for a Serkowski fit. Our new data, which 
have better wavelenth coverage and better S/N, show that AAfa2:~4775 A for S78, very different 
from that found by Martin k. Shawl (1982). Their value of 3500 A is at the extreme low end of 
Xmax values measured in the Galaxy (e.g., Serkowski et al. 1975). The unusual sightline toward 
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HD 210121, which has one of the lowest Ry values ever measured, has Xmax ~3800 A (Larson, 
Whittet, & Hough 1996). None of the ~20 sightlines measured in the LMC show Xmax values 
below 4900 A (Clayton et al. 1983). Only a few sightlines have been observed in the SMC and the 
signal-to- noise is relatively low but Xmax lies in the normal Galactic range (Rodrigues et al. 1997). 

The fits for the four sightlines in M31 reveal values of Xmax ranging from 4800 to 5500 A. 
These are very average values of Xmax as measured in the Galaxy and the Magellanic Clouds. It 
has been found that Ry ~ ^-^Xmax (Serkowski et al. 1975, Whittet &; van Breda 1978; Clayton Sz 
Mathis 1988). The range of Xmax measured in M31 corresponds to Ry values of 2.7 to 3.1. The 
range in Ry seen in the Galaxy is 2.5 to 5.5. So, at least for this very small sample, M31 does not 
show evidence for sightlines with greater than average grain size, if the nature of the grains is the 
same as for the Galaxy. This result is in contrast to the paucity of small grains inferred from IRAS 
colors (Xu & Helou 1994). 

The polarizing efficiency of the M31 dust grains can be estimated using the ratio of peak linear 
polarization to color excess, Pmax/^(B-^)- While the quantity is often equated with the degree of 
grain alignment efficiency, it is also a function of the grain's non-spherical nature as well as the 
angle between the alignment axis and the line of sight. However, for brevity, we simply refer to 
the value as the "efficiency ratio." The maximum value of this ratio found in the Galaxy is 9% 
mag~^ (Serkowski et al. 1975; Clayton & Cardelli 1988). We corrected the estimated reddenings 
for a Galactic foreground of E(B-V)=0.1 mag (Barmby et al. 2000). The values of pmax/E(B-V) 
for our sample are given in Table 1. They range from 7.1 to 15.3% mag~^. These values generally 
support the results of Martin & Shawl (1982) and imply a similar polarization efficiency in M31 and 
the Galaxy. The values of p^aj;/E(B-V) for our sample are all near the maximum value but some 
bias is expected since sightlines known to have significant interstellar polarization were selected for 
inclusion in this sample. The polarization efficiency toward S233 seems very high, but the amount 
of reddening is fairly small and the uncertainty is a considerable fraction of the E(B-V) value itself. 

Baade 327 is of interest as it is both the most luminous globular cluster known and also 
the most reddened globular cluster in M31 (Barmby, Perrett, &: Bridges 2002). van den Bergh 
(1968) suggested that it is strange that the brightest globular cluster in M31 would also be the 
most reddened. This coincidence led to speculation that perhaps the value of Ry might be small 
along this sightline resulting in less total extinction. Ay. Our result for Baade 327, Xmax ~ 5450 
A implies that the dust along this sightline has the Galactic average, Ry ~3.1. Even if Ry = 
2.5, as suggested by Kron &: Mayall (1960), Baade 327 only drops to the second most luminous 
globular cluster (Barmby et al. 2002). They also point out that the most reddened globular 
clusters lie close to the center of M31 and those close to the center of M31 seem to be brighter. 
The polarization efficiency toward Baade 327 is tantalizingly high, particularly given its large color 
excess and its anticipated favorable viewing geometry (i.e., the angle between alignment and line 
of sight) (Martin & Shawl 1979, 1982). Clearly, additional observations of M31 globular clusters 
would provide further insight into the dust environment in that galaxy. 
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Fig. 1. — Spectropolarimetric observations of M31 globular clusters. The upper panels show the 
foreground corrected polarization, the middle panels show position angle and the bottom panels 
show the spectra. Open squares represent the MMT data,the stars represent the WIYN data and 
the triangles represent the Martin k, Shawl (1982) data. The error bars are la. The lines drawn 
represent the best Serkowski fits to the MMT data only. See text. 



Fig. 2. — The Serkowski fit parameters, K and \max ^''^^ plotted for the four M31 globular clusters. 
The error bars are la. The line represents the K = l.QQXmax relationship that holds for Galactic 
interstellar polarization (Whittet et al. 1992). 
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Table 1. Observed M31 Globular Clusters 



Cluster^ 


yb 


E{B - V)^ 


Pmax{%Y 






P^„JE(B-V)(% mag-i)<i 


S78 


14.26 


0.36±0.04 


1.96±0.03 


4775±250 


1.32±0.28 


7.5±0.15 


S150 


15.44 


0.48±0.06 


3.19±0.05 


5200±225 


1.78±0.46 


8.4±0.16 


S233 


15.38 


0.23±0.05 


1.99±0.03 


5275±175 


2.38±0.49 


15.3±0.39 


Baade 327 


16.71 


1.38±0.02 


9.14±0.12 


5450±175 


2.22±0.35 


7.1±0.02 



^Vetesnik (1962); Sargent et al. (1977) 

^Total reddening as estimated by Barmby et al. (2000) 

^Based on MMT data only 

^Galactic foreground reddening of E(B-V)=0.1 subtracted. 



Table 2. S78-MMT 



A-bin (A) 




(% 


) 


1 


9 




Vcorr 


dcorr 


4125 


1.78± 0.36 


41.5 


± 


5.4 


1.86 


35.76 


4375 


2.00 


± 


0.09 


51.8 


± 


1.5 


1.94 


46.32 


4625 


1.89 




0.08 


50.0 


± 


1.1 


1.86 


44.09 


4875 


1.97 




0.05 


49.3 


± 


0.9 


1.95 


43.58 


5125 


2.04 




0.04 


51.5 




0.6 


1.99 


45.89 


5375 


2.01 




0.04 


51.5 




0.6 


1.96 


45.81 


5625 


1.90 




0.04 


52.8 


± 


0.7 


1.83 


46.72 


5875 


1.96 




0.04 


51.5 




0.4 


1.91 


45.63 


6125 


1.92 




0.03 


51.9 


± 


0.5 


1.87 


45.98 


6375 


1.76 


± 


0.03 


51.8 


± 


0.7 


1.71 


45.41 


6625 


1.75 




0.03 


53.6 




0.4 


1.68 


47.25 


6875 


1.77 




0.04 


53.3 




0.5 


1.70 


47.12 


7125 


1.60 




0.04 


52.2 




0.6 


1.55 


45.38 


7375 


1.50 




0.05 


51.9 




0.8 


1.46 


44.79 


7625 


1.54 




0.06 


50.9 


± 


1.0 


1.51 


U.U 


7875 


1.49 


± 


0.05 


52.2 


± 


1.0 


1.44 


45.29 


8125 


1.45 


± 


0.07 


53.0 




1.4 


1.39 


46.09 
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Table 3. S78-WIYN 



A-bin (A) 


p(%) 


e 


Vcorr 


Ocorr 


5500 


2.37 ± 0.23 


61.7 ± 2.8 


2.18 


57.31 


6156 


1.85 ± 0.19 


46.7 ± 2.9 


1.87 


40.54 


6469 


1.81 ± 0.11 


52.1 ± 1.7 


1.76 


45.74 


6781 


1.82 ± 0.09 


52.4 ± 1.4 


1.76 


46.18 


7094 


1.74 ± 0.08 


52.5 ± 1.3 


1.68 


46.09 


7406 


1.81 ± 0.07 


53.6 ± 1.1 


1.74 


47.58 


7719 


1.60 ± 0.08 


57.0 ± 1.4 


1.49 


50.34 


8031 


1.42 ± 0.08 


61.4 ± 1.6 


1.26 


54.30 


8344 


1.44 ± 0.08 


61.0 ± 1.6 


1.29 


54.16 



Table 4. S150 



A-bin (A) p(%) 

Vcorr 9, 



4125 


2.70 


± 


0.44 


36.7 




5.1 


2.82 


33.12 


4375 


2.60 




0.16 


34.3 


± 


1.6 


2.76 


30.60 


4625 


3.11 




0.09 


36.3 




1.0 


3.25 


33.00 


4875 


2.75 




0.09 


36.4 




0.7 


2.89 


32.68 


5125 


3.12 




0.05 


34.9 


± 


0.6 


3.28 


31.63 


5375 


2.90 


± 


0.07 


35.4 


± 


0.6 


3.05 


31.81 


5625 


2.97 


± 


0.06 


35.0 


± 


0.4 


3.13 


31.51 


5875 


3.21 




0.05 


32.7 


± 


0.5 


3.39 


29.63 


6125 


2.81 




0.05 


35.5 


± 


0.4 


2.96 


31.84 


6375 


2.69 




0.05 


35.4 




0.4 


2.84 


31.64 


6625 


2.62 




0.04 


34.7 


± 


0.5 


2.78 


30.92 


6875 


2.66 




0.04 


34.5 




0.6 


2.82 


30.87 


7125 


2.51 




0.05 


35.3 


± 


0.5 


2.66 


31.43 


7375 


2.35 


± 


0.06 


35.8 


± 


0.8 


2.49 


31.76 


7625 


2.34 




0.06 


34.8 




0.9 


2.49 


30.88 


7875 


2.27 




0.08 


36.8 




1.2 


2.39 


32.75 


8125 


2.18 


± 


0.13 


33.1 




1.1 


2.31 


29.18 
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Table 5. S233 



A-bin (A) p(%) 



4125 


1 


.82 


± 


0.42 


53.6 


zb 


5.7 


1 


.74 


47.76 


4375 


2 


.03 


± 


0.18 


44.9 


± 


2.5 


2 


.06 


39.59 


4625 


1, 


.84 


± 


0.12 


44.8 


zb 


1.8 


1 


.88 


38.82 


4875 


1 


.90 




0.11 


43.5 


zb 


1.6 


1 


.96 


37.67 


5125 


1 


.92 




0.08 


41.3 


zb 


1.2 


2, 


.01 


35.64 


5375 


1 


.94 


± 


0.07 


44.0 


zb 


1.4 


1, 


.99 


38.20 


5625 


1, 


.74 


zb 


0.10 


42.6 


zb 


1.1 


1 


.82 


36.21 


5875 


1, 


.78 


zb 


0.07 


44.6 


± 


1.1 


1 


.83 


38.25 


6125 


1, 


.85 


± 


0.05 


43.4 


zb 


0.9 


1 


.91 


37.48 


6375 


1 


.84 




0.06 


43.0 


zb 


1.0 


1, 


.91 


37.11 


6625 


1 


.74 




0.06 


44.5 


zb 


1.0 


1, 


.79 


38.25 


6875 


1, 


.63 


± 


0.07 


47.2 


zb 


1.2 


1 


.64 


40.57 


7125 


1, 


.62 


zb 


0.06 


44.0 


zb 


1.3 


1 


.67 


37.53 


7375 


1, 


.52 


zb 


0.09 


45.5 


± 


1.4 


1 


.56 


38.63 


7625 


1, 


.36 


zb 


0.08 


52.3 


zb 


1.8 


1 


.31 


44.58 


7875 


1 


.34 


zb 


0.14 


37.6 


zb 


2.4 


1, 


.47 


30.82 


8125 


1 


.02 


zb 


0.15 


51.6 


zb 


3.8 


1, 


.00 


41.77 
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Tablc 6. Baado 327 



A-bin (A) p(%) 



4125 


5.72 


± 


3.26 


74.0 


zb 


9.7 


5.40 


73.01 


4375 


6.49 


zb 


1.68 


45.8 


± 


5.2 


6.50 


44.10 


4625 


7.31 


lb 


1.05 


54.5 


zb 


3.3 


7.19 


53.02 


4875 


8.81 


± 


0.74 


58.3 


zb 


1.9 


8.64 


57.10 


5125 


9.09 




0.45 


53.0 


zb 


1.3 


8.99 


51.78 


5375 


9.37 




0.31 


56.1 


zb 


0.8 


9.21 


55.86 


5625 


9.57 




0.23 


53.9 


zb 


0.7 


9.46 


52.73 


5875 


9.18 




0.20 


53.1 


± 


0.6 


9.08 


51.85 


6125 


9.04 


± 


0.18 


53.8 


zb 


0.4 


8.93 


52.58 


6375 


8.63 




0.14 


54.7 


zb 


0.5 


8.51 


53.42 


6625 


8.62 




0.15 


54.9 


zb 


0.4 


8.50 


53.63 


6875 


8.37 


± 


0.13 


55.2 


zb 


0.4 


8.25 


54.00 


7125 


7.74 


zb 


0.13 


54.9 


zb 


0.5 


7.62 


53.59 


7375 


7.33 


zb 


0.13 


53.4 


± 


0.5 


7.23 


52.04 


7625 


7.28 


zb 


0.12 


54.4 


zb 


0.5 


7.17 


53.05 


7875 


7.15 


zb 


0.19 


54.0 


zb 


0.5 


7.05 


52.61 


8125 


6.72 


zb 


0.22 


55.1 


zb 


0.8 


6.60 


53.66 



